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We consider ferroelastic first-order phase transitions with No p order-parameter strains entering Landau free 
energies as invariant polynomials, that have Nv structural-variant Landau minima. The total free energy in- 
cludes (seemingly innocuous) harmonic terms, in the n = 6 — No p non-order-parameter strains. Four 3D tran- 
sitions are considered, tetragonal/orthorhombic, cubic/tetragonal, cubic/trigonal and cubic/orthorhombic unit- 
cell distortions, with respectively, Nop — 1, 2, 3 and 2; and Nv = 2, 3, 4 and 6. Five 2D transitions are also 
considered, as simpler examples. Following Barsch and Krumhansl, we scale the free energy to absorb most 
material-dependent elastic coefficients into an overall prefactor, by scaling in an overall elastic energy density; 
a dimensionless temperature variable; and the spontaneous-strain magnitude at transition A << 1. To leading 
order in A the scaled Landau minima become material-independent, in a kind of 'quasi-universality'. The scaled 
minima in iVop-dimensional order-parameter space, fall at the centre and at the Nv corners, of a transition- 
specific polyhedron inscribed in a sphere, whose radius is unity at transition. The 'polyhedra' for the four 3D 
transitions are respectively, a line, a triangle, a tetrahedron, and a hexagon. We minimize the n terms harmonic 
in the non-order-parameter strains, by substituting solutions of the 'no dislocation' St Venant compatibility 
constraints, and explicitly obtain powerlaw anisotropic, order-parameter interactions, for all transitions. In a 
reduced discrete- variable description, the competing minima of the Landau free energies induce unit-magnitude 
pseudospin vectors, with Nv + 1 values, pointing to the polyhedra corners and the (zero- value) center. The total 
scaled free energies then become Zm v +i clock-like pseudospin hamiltonians, with temperature-dependent local 
Landau terms, nearest-neighbor Ginzburg couplings, and powerlaw St Venant interactions that drive the elastic 
domain-wall texturing. The scaled free energies can be used in relaxational or underdamped dynamic simula- 
tions, to study ferroelastic strain textures and their dynamical evolution pathways. The pseudospin models can 
similarly be studied via local meanfield treatments, and Monte Carlo simulations. 



I. INTRODUCTION 

Although first order phase transitions predominate in na- 
ture, second-order transitions have attracted much theoretical 
attention, because of the taxonomically simplifying concept of 
'universality classes'—. Diverse materials in the same univer- 
sality class have common critical exponent values as a scaled 
temperature T/T c — > 1, that depend only on the dimensional- 
ity d of coordinate space; the number Nop of order parameter 
components; and the number Nv of degenerate energy min- 
ima or 'variants'. Spin models are prototypical, with d = 2, 
Nop = 1, Nv = 2 for a 2D Ising model. The symmetry- 
breaking transition is signalled by spontaneous nonzero aver- 
ages of a spin component, as one of the degenerate-minima 
variants in order parameter {OP) space is picked out. First- 
order transitions by contrast, seem to be inherently material 
specific, as they lack a divergent length scale to induce univer- 
sality, by rendering irrelevant the finite-scale material parame- 
ters. On the other hand, precisely because critical fluctuations 
are unimportant, an approach focusing on free energy minima 
is more reliable. If the free energies could be scaled to make 
at least these minima independent of the material coefficients, 
one would have a kind of 'quasi-universality' for first-order 
transitions. 

Ferroelastic, displacive structural transformations 2 ^ as 
in the austenite/martensite transition are (mostly) first or- 
der, symmetry lowering transitions, with the discrete sym- 
metries of high/low temperature lattices often having a 
group/subgroup relationship. Lattices in d dimensions have 



d(d + l)/2 Cartesian strain-tensor components, whose lin- 
ear combinations are the same number of physical strains, 
of which Nop are the order parameters. The Landau vari- 
ational free energies are sums of high-order polynomial in- 
variants in the OP, with many material-specific anharmonic 
elastic coefficients^. 7 , that take effort to extract from experi- 
ments. With a single high temperature zero-strain state, and 
Nv other structural variants apppearing as temperature is low- 
ered, there are Ny + 1 degenerate Landau minima at the first- 
order transition temperature. Twinning, or spatial coexistence 
of competing structures separated by oriented domain walls, 
is common 2,8 . The orientation comes from long-range elastic 
forces^ or powerlaw anisotropic interactions 11-17 between the 
No p order parameters, that are present but hidden, in the con- 
ventional displacement representation. As will be seen, these 
effective interaction*^, arise from a constrained minimization 
of free energy harmonic terms in the n = d(d + l)/2 — Nop 
non-order-parameter strains, subject to St. Venant compatibil- 
ity conditions'^. 

Barsch and Krumhansl (BK)i£ have scaled the Landau free 
energy of a 2D square/rectangle transition, and the 3D cu- 
bic/tetragonal transition using three scaling parameters, to a 
dimensionless form that is internally independent of elastic 
coefficients. A conceptually important dimensionless temper- 
ature variable r(T) = (T - T C )/(T - T c ) can be introduced, 
that absorbs elastic coefficients of the quadratic term in the 
OP-strain magnitude. It is unity t(To) = 1 at the first order 
transition temperature T = Tq; while it vanishes r(T c ) = at 
the lower spinodal T c , 
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In this paper, we generalize this BK procedure for the 
Landau free energies of four 3D and five 2D transitions, 
that can have Nop — 1,2,3 order parameters, and Ny = 
2,3,4,6 variants, to absorb the (often unknown) material 
coefficients^' 7 , into an overall energy prefactor. The four 3D 
transitions are tetragonal to orthorhombic, cubic to tetrago- 
nal, cubic to trigonal; and cubic to orthorhombic, with their 
number of unsealed material coefficients respectively N mat = 
3,3,4 and 6. The 3D transitions we have chosen are rele- 
vant for functional materials : tetragonal/ orthorhombic (as 
for high-temperature superconductors such as YBCO); cu- 
bic/ tetragonal (as for shape memory alloys such as FePd); 
cubic/orthorhombic (as for ferroelectrics such as BaTiO); 
and cubic/trigonal (as for CMR oxides such as LaSrMnO). 
We show that, apart from the overall energy-density coeffi- 
cient, the scaled free energies evaluated at the Ny + 1 min- 
ima are material-independent in the first three 3D cases, and 
weakly dependent on a single material parameter in the cubic- 
orthorhombic case. The minima fall at the center and corners 
of a transition-specific polyhedron inscribed in a sphere, of 
unit radius at transition. The 'polyhedron' for Nop = 1,2 
or 3 can be a line, a triangle, a hexagon or a tetrahedron. We 
evaluate the St Venant compatibility potentials through their 
Nop x Nop matrix kernels, for all the nine transitions us- 
ing a constraint substitution method, that is more direct than 
a Lagrange multiplier 11 method. The scaled total free ener- 
gies can be used in over-damped, or under-damped OP strain 
dynamics (that includes Langevin noise terms with powerlaw 
anisotropic spatial correlations)i2iI2. In a reduced, discrete- 
variable description that retains only the Landau minima, the 
total scaled free energies induce clock-!£ like 'Zjv v +i' hamil- 
tonians in terms of unit-magnitude, Nop -dimensional vec- 
tors, pointing to Ny + 1 values. The hamiltonians are bilinear 
in the pseudospins, with (temperature-dependent) quadratic 
on-site contributions from the Landau term; nearest-neighbor 
ferromagnetic couplings from the Ginzburg term; and power- 
law anisotropic interaction potentials from the St Venant term. 
The pseudospin hamiltonians can be used for local meanfield 
treatments, and Monte Carlo simulations, of ferroelastic tex- 
tures. The models are also be relevant for complex functional 
materials, with lattice strains coupled to intracell charge, spin 
and orbitalsii 

In more detail, the transition-specific scaling procedure in- 
volves a choice of three scaling parameters: a typical spon- 
taneous strain magnitude A that is small; a typical elastic en- 
ergy density Eq\ and a (Landau) first-order transition temper- 
ature To chosen such that the scaled BK temperature variable 
t(Tq) = 1, and the scaled OP magnitude is unity at transition. 
The smallness of A << 1 justifies a finite-sum truncation of 
the expansion in invariants; and a neglect of geometric non- 
linearities inside the Lagrangian strains, as a perturbative first 
approximation. Working to leading order in A, yields a 'quasi- 
universal' scaled Landau term that has an overall, material- 
dependent energy density Eq, as mentioned. The total free 
energy has Ginzburg and St Venant terms that determine elas- 
tic domain-wall texturing, and carry (two) material-dependent 
coefficients. Elastic pseudospin models in the context of 
martensites have been considered by several groups^—, but 



the models obtained here are explicitly induced by scaled free 
energies. 

In Sec. II we set up the general BK scaling procedure 
for 2D and 3D free energies. We then consider transitions 
in increasing number of variants, Ny = 2,3,4,6. In Sec. 
Ill we consider two-variant, single order parameter Ny — 
2, Nop = 1 cases, namely the 3D tetragonal/orthorhombic 
transition; and in 2D the square/rectangle or square/ rhombus; 
and rectangle/oblique-polygon transitions. Section IV con- 
siders cases with Ny = 3 and Nop = 2, namely the 3D 
cubic/ tetragonal, and 2D triangle/ centered- rectangle transi- 
tions. Section V considers Ny = 4 cases: the 3D cubic/ trig- 
onal (Nop = 3), and 2D square/ oblique (Nop = 2). Section 
VI presents the Ny = 6, Nop = 2 case of 3D cubic/ or- 
thorhombic and 2D triangle/ oblique transitions. Section VII 
obtains the reduced pseudospin hamiltonians, for all the tran- 
sitions. Section VIII outlines possible simulation approaches, 
while the last Section IX is a summary. An overview in Table 
I collects the generic numbers of the scaled Landau free en- 
ergies, common to all materials with a given transition in the 
same quasi-universality class. The Appendix derives compat- 
ibility kernels for the four 3D and five 2D transitions consid- 
ered, through the direct substitution method, contributing to a 
'library of kernels' for use in simulations. 

II. SCALING PROCEDURE 

Here (A) we define elastic variables, and (B) state the gen- 
eral scaling procedure. 

A. Distortion and strains 

The distortion tensor- £> M „ can be defined in terms of gradi- 
ents of the displacement vector u{r) of points in a continuum 
medium, by 

D liv {f)=du^(f)/dr u , (2.1a) 
where fx, v run over x, y, z; or in Fourier space, 

D liV (k)=iu lt (k)k v . (2.16) 

The D tensor is a sum D = e + w of a symmetrized (e) 
and antisymmetrized (w) distortion (or local rotation) tensor. 
With "T" a transpose, 

e(f} = i[D + D T ]; w(f) = \[T> D T ]. (2.2) 

We will refer to the symmetrized distortion tensor e^„ as the 
'strain' tensor, and take it as the working variable, in a strain 
representation. It is distinct from the 'Lagrangian-strain' E^, 
that is a derived quantity* defined below. 

In the nonuniform case, the six symmetrized distortions 
e All/ (r) cannot vary arbitrarily, if lattice integrity is to be main- 
tained (i.e. no defects such as dislocations). They are linked 
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by the St. Venant's 'compatibility' equations 9,10 of 1864 that 
ensure distorted unit cells fit together in a smoothly compati- 
ble way. With a double-Curl operation defining an Incompat- 
ibility or 'Inc' operation, in coordinate and Fourier space, the 
Cartesian-component constraints on e(f) are 



Inc[e(r)} = V x (V x e(r)) T = 0; (2.3a) 



Inc[e(k)} = kx e{k) xk = 0. (2.36) 

(Similar equations hold for w(r), but for the transitions con- 
sidered here we take spontaneous local rotations to be zero.) 

From (2.1), in the displacement representation where u(f) 
is the working variable, compatibility is satisfied as an identity 
since Curl(Grad) = 0. Baus and LovettiS have proposed the 
distortion tensor D be taken as the working variable, with the 
St Venant compatibility equations then being field equation 
constraints to be satisfied, rather than identities. In an electro- 
magnetic analogy 10 this is like working with the magnetic in- 
duction B rather than the vector potential A, with DivB = 
then a Maxwell field equation to be satisfied, rather than a 
Div(Curl) = identity. The change to a distortion-tensor 
variable is natural and useful, since the free energy depends 
directly on the distortion; and the constrained minimization 
in the distortion reveals powerlaw anisotropic potentials^— , 
that are hidden in the more conventional displacement repre- 
sentation. 

The elastic free energy is invariant under global uniform 
rotations 5 . Consider a line-element in the elastic medium, de- 
scribed by a small imbedded column-vector a at a site r. Un- 
der distortion of the medium, it locally stretches or rotates to a 
vector A(r) = [1 + D(r)]a. Free energies F can depend only 



on the scalar product A (r) .A(r) 



[1 + D] y [l+D] 



a T .[l+2E].a, where A T , a T are row vectors. The free energy 
F{E I1V ) thus depends on components of the 'Lagrangian- 
strain' tensor— E^, that is a derived quantity, defined in terms 
of the basic distortion variable by 



E(r)^l[D 5 



D] + Jd t D 



e(f)+g(r), (2.4) 



tD D. 



where the 'geometric nonlinearity' is g _ 

It is convenient to consider d(d + l)/2 physical strains e a 
that are linear combinations of the d(d+ 1)/2 Cartesian strains 
e^ v . For 2D lattices, with three Cartesian components, they 
are e a with a = 1,2,3: the compressional (ei), deviatoric 
(e2), and shear (e3) strains. With X, Y transforming as 2D 
Cartesian coordinates, the physical strains transform as e\ ~ 
X 2 + Y 2 , e 2 - X 2 - F 2 , e 3 - XY, and are defined asi^ 



The physical Lagrangian-strains are similarly defined, so 

Ei{r) = -j( E xx + E yy ) = e x + g x \ 



E 2 (r) = ~^{E XX - E yy ) = e 2 + g 2 ; 



E 3 (f) = c 3 E xy = e 3 + 53; (2.56) 

where the physical geometric-nonlinearities are gi = 
t(9xx + St/7/), 02 = y(9xx ~ 07,7,), 03 = c 3 g xy . For square 
unit-cells in the high-temperature phase, the normalizing co- 
efficients are chosen as c\ = v2 = c 2 , and c 3 = 1, while for 
triangular unit cells C\ = c% = c 3 = 1, 

For 3D cubic lattices with six Cartesian strain compo- 
nents, the physical strains e a with a = 1,2, ..6 are the 
compressional (ei), deviatoric (e 2 , e 3 ), and shear (e$, e^, e^) 
symmetrized-distortions. These transform as one, two, and 
three dimensional irreducible representations of the cubic 
point group. For the cubic lattice, with X, Y, Z transform- 
ing as 3D Cartesian coordinates, the physical strains trans- 
form as ei ~ X 2 + Y 2 + Z 2 ,e 2 ~ X 2 - Y 2 ,e 3 ~ 
X 2 + Y 2 - 2Z 2 ,e 4 ~ YZ,e 5 ~ ZX, e 6 - XY, and are 
defined asiii£ 



ei(r) = ^( e xx + e yy + e zz ); 



e 4 (r) = 2e yz , e 5 (r) = 2e zx , e 6 (r) = 2e xy . (2.6a) 

The physical Lagrangian-strains E a are similarly defined so 
e.g. Ei = (E xx + Eyy + E zz )/V3 = e 1 + g x where g 1 = 

(gxx + 07,7, + gzz)l \/3, and so on. 
For the tetragonal lattice 16 , 



ei(r) = -j={e xx + e yy )- 



1 



62(f) — ~^(g xx 



e 3 (r) 



(2.66) 



ei(f) 



ci 



C2 



(e X x + e yv ); e 2 (r) = — (e 



--yyh 



and the physical shears are the same as (2.6a). These trans- 
form as irreducible representations of the tetragonal point 
group. 

As noted in the Appendix, the St Venant constraint (2.3b) 
in Fourier space can be written in terms of the physical strains 

as 



C3 



e 3(r) = TT^xy + C-yx) = c 3 e xy . (2.5a) 



E O^e a (k) = 0, (2.7) 



a=l,2...6 
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where Oa \k) are compatibility coefficients appropriate to 
the symmetry, labelled by shear components, s = 3 in 2D (one 
constraint), and s = 4, 5, 6 in 3D (three constraints). This is 
used in the Appendix to calculate the compatibility potentials. 

The physical strains can be separated into No p order pa- 
rameter strains labelled by a = i or {e^}; and n = [\d{d + 
1) — Nop] non-OP strains, labelled by a = i or {e^}. (An 
associated separation of the physical Lagrangian-strains for 
each transition into {£/} and {Ei} is also made.) For a fer- 
roelastic transition, the free energy F = Fl + F non + Fg 
. The Landau part of the free energy Fl(E^) can be writ- 
ten as a sum of p-th order polynomial invariants in the OP 
Lagrangian-strains, that can be found by direct evaluation^— 
for simple cases. The non-OP free energy F non (Ej) is taken 
as harmonic in the non-OP Lagrangian-strains. For simplicity 
we neglect symmetry-allowed anisotropic gradient s 14 ' 15 , and 
consider only OP gradient-squared costs in the Ginzburg term 
FoiyEi). For all 2D transitions, an exhaustive evaluation of 
the OP and OP-gradient invariants, and allowed OP/non-OP 
couplings, was carried out through the program ISOTROPY 
of Stokes and Hatchi^. 

Then with a sum or integral (J2r I d d r/ao) over all 
positions f where a is a lattice scale, the total variational free 
energy in terms of free energy densities / is 

F = 2 h ^ + /g(V^) + f non {Ei). (2.8) 

T 

The Landau term Jl is not just an arbitrary Taylor series, but 
is a finite-sum expansion in terms of symmetry-allowed, in- 
variant polynomials of physical OP Lagrangian-strains 

f L = C^(T)I 2 + °pC {p) I v {E t ). (2.9a) 

p=3,...p max 

The second order invariant I 2 = £\ Ef is common to all 
transitions, and is separated out. The anharmonic elastic co- 
efficients > are temperature-independent, and a p = 
+1, — 1 are chosen signs, to get Ny minima. We consider an 
elastic coefficient C { ^(T) = (T - T C )C (2) that would par- 
tially soften to zero on cooling to a temperature T = T c , but 
is preempted by the first order transition. The Ginzburg and 
non-OP terms, with domain wall cost parameter b and A^ l > 
elastic coefficients^, are 

f G = Y,KVE e (r)) 2 ; fnon = £ l^Efif). (2.9b) 
i i 1 

The temperature-dependence of sound velocities 2 - 7 - experi- 
mentally determine the linear slope C^ 2 \ and the extrapo- 
lated temperature T c . The curvature of phonon spectrum at 
long wavelengths 27 determine b, that is related to domain-wall 
energy costs at short wavelengths. We take it to be positive 
b > 0; but for b < 0, one can add phenomenologically a 
symmetry allowed positive-coefficient, fourth-order gradient 
term, for stability. 



B. General scaling and minimization procedure 

We now follow a Barsch-Rrumhansl procedure^: scal- 
ing all Cartesian distortions £)^„ — > XD^ or — > 
Xe^y , — > XWf^u in a typical value A; scaling all energy 
terms in Eq ; and defining a transition temperature To- These 
three parameters are chosen in terms of the material-specific 
elastic constants, to make the Landau free energy simple. The 
parameters have physical meanings: since scaled terms are of 
order unity, the overall prefactor Eq is the elastic energy per 
unit cell; the temperature Tq(> T c ) is the first-order transition 
temperature that pre-empts the second-order elastic-constant 
softening at T c ; and A is the spontaneous-strain magnitude at 
To, since the scaled strain is chosen to be unity at transition. 

Since the physical distortions are linear combinations of the 
Cartesian distortions, they change as e a — > Xe a , and the phys- 
ical Lagrangian-strains change as 

E a -> XE Q {X) = A(e Q + Xg a ) . (2.10) 

Henceforth e a is the scaled symmetrized-distortion or scaled 
strain. The (scaled) geometric nonlinearity g a carries a pref- 
actor A. The free energy changes as 

F{E a ) -> F[XE a {X)] = E F[E a (X)], (2.11a) 

defining a dimensionless F[E a (X)] — J2f fi^a(X)] where 
the scaled dimensionless free-energy densities are 

f[E a (A)] = /l[^(A)]+/ g [V^(A)]+/„ „[^(A)]. (2.116) 

These terms contain dimensionless scaled coefficients that ab- 
sorb the 'external' A powers of (2.10) as 

C p = Cto\*/Eo; Ai = A^X 2 /E ; £ 2 = &A 2 /ao 2 £ . (2.12) 

Here a dimensionless length £ has a lattice constant ao scale, 
from a substitution as below in the Ginzburg term V — > A/ao, 
where A is a discrete-difference operator on a computational 
grid. 

A scaled temperature^ r absorbs the harmonic-term mate- 
rial dependence, and is defined as 

r(T) = (T-T c )Co i2] X 2 /E = (T-T C )/{T -T C ), (2.13) 

where the transition temperature T > T c is determined by 
requiring r(T ) = (T - T C )C (2) X 2 /E = 1 or 

To = T c + S /(C (2) A 2 ). (2.14) 

At the first-order Landau transition temperature at a uni- 
versal value t(Xo) = 1, the nontrivial or 'martensite' wells 
are degenerate with the trivial or 'austenite' well; while T c is 
the lower spinodal, where at a universal value r(T c ) = 0, the 
metastable austenite well disappears. The upper spinodal T up 
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where the martensite-variant wells disappear, turns out also to 
have a universal value r{T up ) = r up for three of the four 3D 
transitions. 

A constrained minimization of the scaled free energy 
f[E a (X)] of (2.11b) with respect to e a , would, for general 
A, fix the parameters A, Eo, and yield an effective OP-OP in- 
teraction potential. However because of geometric nonlineari- 
ties the calculation is involved, and we perturbatively evaluate 
the scaling parameters and compatibility potentials, as below. 
Corrections can in principle, be calculated. 

We make a simplifying assumption that the typical sponta- 
neous distortion is small compared to unity 

A«l, (2.15) 

and for most materials this is indeed a few per cent, A ~ 10~ 2 . 
Then the scaled physical Lagrangian-strains E a (X) = e a + 
Xg a are approximated by the scaled physical symmetrized- 
distortions, 

E e {X) ~ E t (0) = e r ,Ei(X) ~ ^(0) = e*. (2.16) 

In fact, this approximation of dropping geometric nonlineari- 
ties in the Lagrangian-strains is commonly made without spe- 
cific comment. In the conventional (unsealed) displacement 
representation it is implicitly justified as a long-wavelength 
truncation: in Fourier space the strain tensor is ~ [cfu(g)] M „, 
while the geometric nonlinearity is g^iq) ~ [qu(q)] 2 , that 
is higher order in q — > 0. Instead, in the (scaled) strain repre- 
sentation, the neglect of geometric nonlinearities is seen as the 
leading term in a small-parameter expansion in A, that could 
be systematically corrected. 

The scaled free energy densities of (2.11b) then become 
f[E a (X)] ~ f{e a ), to leading order in A, separately for each 
distinct symmetry-invariant term, 

f(e a ) = f L (e e ) + f G (Ae e ) + f non (eA, (2.17) 
From (2.9a), (2.9b), (2.12) and (2.13), 

f L (e i ) = (T-l)^e i 2 +f {e i ), (2.18a) 

t 

/o(e £ ) = ^e c 2 + ]ra p C p / P M, (2.186) 

i V 

/ G (Ve,H^ 2 (Ae,) 2 , /„ on ( e «) = EY e -( 2 - 18c ) 

e i 

Here as mentioned previously, A = (A^, A y , A z ) has dis- 
crete forward-difference operator components on a cubic 
computational grid. 

We pause to relate the unsealed harmonic coefficients A^ % > 
of (2.9) to the material elastic constants. The elastic en- 
ergy is* 6 - E = j J2 a 8 C a px a Xj3 where with a = 1, 2, ..6, 



the {x a } are Cartesian strains written as a column vector 
such as xi = e xx ,X4 = 2e yz . In the cubic case, there 
are three independent elastic constants in the Voigt nota- 
tion, Cn, C12, C$ji{= C55 = 6*66)- Writing Cartesian in 
terms of physical strains, through (2.6a) and (A9), the en- 
ergy E is diagonalized, and a comparison with (2.9) yields 
C n - C12 = C^(T),Cn + 2C 12 = A«,C 4 4 = A^\ 
It is useful to introduce the elastic anisotropy parameter 26 
A(T) = 2C Ai /{C lx - C12), where A > 1 (or A < 1) cor- 
responds to greater stiffness in the body diagonal < 111 > 
directions (or cubic axis < 100 > directions). Strongly 
anisotropic materials can have A ~ 10. Then from (2.12), 
the scaled shear coefficient Aj(= As = Aq), and the scaled 
compression coefficient A\ are both in terms of the elastic 
anisotropy parameter A(T = To) at transition, 

A 4 = A(T ); A 1 = jA{T ); 7 = [{Cn+2C 12 )/C M ]. (2.18d) 

The elastic constant ratio 7 = j4i/Ai(= A1/A5 = Ai/Ae), 
that enters the compatibility potentials, can for simplicity 
be set in simulations to a constant, say A^ /2A^ 4 > ~ 1 
as inii FePd. [For xy plane distortions, the shear term is 
(Ae/2)4e xy 2 , so from (2.5a) the scaled 2D shear coefficient 
appearing later is also proportional to the anisotropy at tran- 
sition.] Finally we note that from (2.12), the scaled Ginzburg 
coefficient is similarly £ 2 ~ (b/al)A(T a ). 

We now return to the main argument. A constrained min- 
imization of the harmonic non-OP terms in (2.18c) as in 
the Appendix yields the non-OP in terms of the OP strains, 
e i(^) = Bu(k)ei(k) where the coefficients are in 
terms of the coefficients a ^ 8 ' of (2.7). Substituting back in 
the harmonic term, F non [ei(ee)} = F compat {ee) induces the 
St Venant term 

F comp at(ee) = ]T U U '(k)e e {k)e* e ,{k). (2.19) 

k,£,t' 

The compatibility kernel AiUu'(k) = Y^iAiBig(k)B*g, (k) 
is evaluated for each transition in the Appendix, and is essen- 
tially dependent only on the wave-vector direction fc, indepen- 
dent of the magnitude \k\. In coordinate space, the compati- 
bility potential is hence an anisotropic powerlaw with a falloff 
exponent equal to the dimensionalit y 1 1 ~ 17 - 23 JJ U , (i?) ~ 1/ R d . 
(Write the Fourier integral of Uw (fc)and change the wave- 
vector integration variable |fc| — > \k\/R: the exponent simply 
comes from the phase space dimension.) 

We focus on the Landau term Jl. To find the minima in OP 
space it is convenient to work in polar coordinates, following 
Toledano and Toledano 6 -. For example the Nop = 2 dimen- 
sional 'vector' in OP space is e = (e cos <f>, e sin cf>), where the 
'radial' variable is 

*s|S| = [5> a ] 1/a , (2.20) 

t 

and the Landau free energy density is 



6 



/i(e) = f L (e, J>) = (r - l)e 2 + / (e, 0), (2.21) 

where the transition-specific /o is temperature independent. 
We demand that the nontrivial Landau minima are at m = 
l,2,..Ny equivalent points {e m ,0 m }, with the same radii 
e m = e in OP space. The conditions are: 



= 2(r-l)e f 



dfo{e m , 4>m) 
de m 



= 0, (2.22a) 



locating the martensitic minima in the radial direction at a 
temperature-dependent e = e m = £{t); and 



df L (e r 



I 



d<t>„ 



0, (2.226) 



locating the minima in the azimuthal direction at a 
temperature-independent <p = 4>m- At transition r = 1, we 
also demand that the nonzero minima on a Nop -dimensional 
'sphere' of radius e(t = 1) = 1, become degenerate with the 
trivial minimum fp{e = 0) = 0. Hence 



fo(e 



1,< 



0. (2.22c) 



Above an upper spinodal t = r up the radial solutions e(r) 
become imaginary, and there is only the trivial austenite min- 
imum. It is convenient for later use to define the Landau free 
energy at minima 

f L ( T ) = f L (e, cj> m ) = £(t) 2 9l (t), (2.22d) 

where gi (r) changes sign at the Landau transition. 

With r defined, we choose the remaining two scaling pa- 
rameters A, Eq so the fo conditions of (2.22b), (2.22c) are 
satisfied. It is useful from (2.21) to separate the angular de- 
pendence into a part A/ (e, (f>) = fo{e, <j)) — fo(e, <f) m ) that 
vanishes in the minimum angular directions, so 

f L (e, 0) = [(r - l)e 2 + f (e, J> m )] + A/ (e, 0), (2.23a) 

and e(r) is determined through minimization of only 

/l(£, <t> m ) = [(t - l)e 2 + /o(e, 4> m )}- (2.236) 

There is always an overall material constant Eq for the 
Landau energy fa = Eq/l, that absorbs unknown higher- 
order elastic coefficients, and can be treated as a fitting pa- 
rameter. Since we work only to leading order in A, any 
material-independence found in the scaled /^(e) contribu- 
tions is strictly speaking only quasi-universal. Landau quasi- 
universality can be of three kinds: (i) strong, i.e. the scaled 
/z,(e, <f>) is independent of material parameters for all e, <f>; 
(ii) medium, i.e. material coefficients appear only in A/o 
that vanishes at <f) = 4> m , so along minima angles, 4> m ) 



is material-independent; and finally (iii) weak, with residual 
material-dependence in fi,{e, 4> m ), and hence in the OP mag- 
nitude e(t) (that is however still unity at transition for all 
materials). The tetragonal/orthorhombic and cubic/ tetrag- 
onal; the cubic/trigonal ; and the cubic/orthorhombic tran- 
sitions (with unsealed material coefficients of respectively 
Nmat = 3,3,4 and 6), turn out to have quasi-universality 
in Jl of respectively the first, second and third kinds. Table 
I summarizes the generic numbers for all transitions consid- 
ered, with different materials with the same transition falling 
into the same 'quasi-universality class'. 

Going back to unsealed variables denoted by primes, the 
unsealed entropy-density difference relative to the austenite 
from the Landau term is s^'(T) = —df' L (e'(T),<f> m )/dT, 
Since the derivative of (2.22a) with respect to e'(r) vanishes, 
only the explicit r-dependence of fi contributes. The scaled 
entropy-density difference is 

s l (t) = [(To - T c )/E ]s l '(t) = -e(r)\ (2.24) 

and is (minus) unity at transition. Of course there are other 
free energy terms, and hysteresis from domain-wall textures, 
so this is just a formal result. We will consider proper fer- 
roelastic transitions with free energy nonlinearities in the 
OP strain driving the transition (without intracell shuffles); 
with high/low temperature unit-cell symmetries having a 
group/subgroup relationship; and without coupling to other 
fields. (There are also improper ferroelastics, with only har- 
monic terms in strains, that are however coupled to other fields 
such as electric polarization or magnetization, whose nonlin- 
earities can induce a structural transition^.) 

In Sees. III-V we consider scaling of the Landau free en- 
ergy (and other terms), for four 3D transitions with No p = 
1,2,3 and for five 2D transitions with Nop = 1,2, presented 
in increasing number of variants Ny = 2, 3, 4, 6. Some cases 
had been scaled earlie r 12 ' 13 ' 18 but are summarized here for 
completeness. The final results for the minima are summa- 
rized in Table I. 



m. TRANSITIONS WITH N v = 2 

We consider one 3D transition and two 2D transitions, all 
with Ny = 2 low temperature variants and a single or- 
der parameter (OP) component Nop = 1, with the number 
of non-OP variables n = ^d(d + 1) — Nop- The transi- 
tions are: (a) tetragonal/orthorhombic in 3D (n = 5); and 
(b) square/rectangle (that includes square/ rhombus); and (c) 
rectangle-oblique cases (all n = 2). See Figs 1 and 2. 



A. Tetragonal/orthorhombic case in 3D: 

N v = 2, Nop = l,n = 5 

There are two deviatoric distortions &2 , 63 in 3D, and the 
single (No p — 1) order parameter is e2 ~ X 2 — Y 2 that can 
change a tetragonal square-cross-section to an orthorhombic 
rectangular-cross-section. Since there are two possible such 
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rectangular elongations (along mutually perpendicular direc- 
tions), one expects two possible variants, as in Fig 2. The 
tetragonal point group G =PA/mmm with vq = 16 elements 
goes in a symmetry-lowering transition, to the orthorhombic 
subgroup g =Pmmm with v a = 8 elements- . The orthorhom- 
bic group describes symmetries of a unit cell with a particular 
rectangular orientation, corresponding to one variant. The ra- 
tio of the number of elements in a point group to that in a 
subgroup, or dimension of the coset G/g, is an integer, that 
we assume^ corresponds to the number of variants. Thus here 
there are iVy — vc/v g = 16/8 = 2 variants, as expected. 
The ratio of the number of (rotational) elements in the point 
group elements, has been taken to be the number of variants 2 , 
yielding the same result 8/4 = 2 here, and in other cases. (See 
however, Section VI below.) 

For these two variants in a first order transition, we need 
up to p m ax — 6 even-order strain invariants I p . The unsealed 
free energy, with sign choices 04 = — 1,(Tq = +1 is fp = 
C (2) I 2 - C {4) h + C (6) / 6 , where I p = e 2 p are invariants, so 
there are N mat = 3 material constants. With e 2 — > Xe 2 and 
C p = \p_C^/E as in (2.12), the scaled Landau free energy 
density f L = f L /E is as in (2.21): 

/L(e a ) = (r-l)el + /o(e a ), (3.1) 
where the temperature-independent /o is 

/o(e 2 )=el-C 4 ef + C 6 e%. (3.2) 

The conditions (2.22a), (2.22c) for degenerate minima are 
df /de + 2(r - l)e = 0, and / (1) = 0. This fixes the two 
coefficients 

C 4 = 2; C 6 = 1, (3.3) 
achieved by choosing the two scaling parameters as 

A = (C^/2C^ 2 ; E = C^(C^/2C^) 3 . (3.4) 
Then /o becomes a perfect square, and 

h(e 2 ) = {r-l)el + ei{ei~l) 2 , (3.5) 

manifestly showing the triple-minima degeneracy at fp = 
for r = 1. The variant minima are at ±£(t) where the order 
parameter magnitude 

e(r) = ||[l + Vl-3r/4]} ' , (3.6) 

is unity at transition. Here from (2.22d) the variable used later 

in pseudospin hamiltonians is fp/e 2 = gp = t — 1 + (e 2 — 

r ' 1/2 

l) 2 . Barriers at e 6 (r) = i |[1 - y/l — 3r/4] > exist in 

the range < r < r up = 4/3. The barriers merge with the 



metastable martensite (or metastable austenite) minimum at 
the upper spinodal r = 4/3 (or lower spinodal r = 0). 

Since fp [S(t)] is independent of material parameters, this 
is quasi-universality of the first kind, with fp = E$fp having 
only an overall material dependence through Eq, that absorbs 
the higher-order elastic constants. 

With the OP sign formally written as an angle, e 2 = 
\e\ cos <fi where 1 62 1 = |e|, and the minima are at <j> m = 
2(m - l)ir/N v , with m = 1, 2(= N v ), where sin2</> m = 0. 
At transition, the nontrivial Landau minima fall at N\- = 2 
points at ±1 on a line in the Nop = 1 dimensional order pa- 
rameter space, as in Fig 3 and Table I. The number of distinct 
martensite/martensite domain walls between variant pairs is 
N w = 1. 

Had we included an eighth order invariant ~ C^ 8 'e 2 8 
in (3.2), the minimum condition would become r — 2e 2 + 
3e 4 + 4C*8e 6 = 0. Here the new scaled coefficient is Cs ~ 
C^X 8 /E = {C^/C^)\ 2 . Assuming the ratio of eighth 
and sixth order unsealed constants is not too large, the shift in 
the roots arising from Cs ~ A 2 <C 1 is negligible. The ex- 
tra eighth order invariant is thus 'irrelevant', in the sense that 
the polyhedral minima remain essentially unchanged, justify- 
ing the finite-sum restriction to p max = 6 of the polynomial 
expansion. 

The scaled Ginzburg term is fa = £ 2 (Ae2) 2 . There 
are n = 5 non-OP strains, namely the compression ei ~ 
X 2 + Y 2 ; the other deviatoric strain e 3 — X 2 + Y 2 -2Z 2 ; and 
the three shears e±, es, eg. Using the three compatibility con- 
straints of (2.7) to eliminate the shears, and minimizing in e±, 
the non-OP strains are determined by the OP. Substituting into 
the harmonic non-OP terms yields an OP compatibility poten- 
tial term as in (2.19): j n on{e-i) = Efci,3,4,5,6^/ 2 ) e i 2 ~^ 
fcom P at{e 2 ) where 

fcom P at = ^U(k)\e 2 (k)\ 2 , (3.7) 

and the kernel U (k) is given in (A26) of the Appendix. 

B. Square/rectangle case: Nv = 2, Nop = 1, n = 2 

The single 2D deviatoric distortion e 2 ~ X 2 — Y 2 turns 
a square to a rectangle. Since the rectangular elongation can 
be along two axes, one expects two variants, as in Fig 1 . The 
point group G = pAmra for a square unit cell has vq = 4 
elements, while the subgroup g for the rectangle is p2mm 
with v g = 2 elements. Thus the number of variants is 15 
Nv = vc/vg = 2 as expected, and we again need up to 
Pmax = 6 order invariants. The scaled Landau free energy 
fL(e 2 ) is the same form as (3.5) above, so the minima are at 
the same ±e(r) of (3.6). 

The scaled Ginzburg term is fa = £ 2 (Ae2) 2 . There are 
n = 2 non-OP strains e\ ~ X 2 + Y 2 , ~ XY, and the har- 
monic free energy term is / no „(e J ) = \ X)f=i )3 A\ei{k)\ 2 . 

This induces a compatibility kerneliiri^i^ U (k) as in (A4) of 
the Appendix. 
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The square/rhombus transition has Ny = 2, Nop = 
l,n = 2, and the 2D shear physical distortion e 3 ~ XY 
as the single order parameter. However this is not an inde- 
pendent transition 15 , since e 2 and e% interconvert through a 
global rotation of Cartesian axes by ir/4. Nonetheless, as an 
exercise the symmetry group of the square is G = pAmm with 
vq = 4 components, while the rhombus symmetry is g = P2 
with v g = 2 elements, so Ny = 2. The scaled Landau free 
energy is Jl{&z), formally as in (3.5), 

f L (e 3 ) = (T-l)ej + el(el-lf,(3.8) 

with two nonzero minima at the same values es = ±S(t). 

The Ginzburg term is f G = £ 2 (Ae 3 ) 2 . The n = 2 
non-OP strains are now compressional and deviatoric, e\ ~ 
X 2 + Y 2 and e 2 ~ X 2 - Y 2 , with the harmoniciii 3 . f non = 
i J2i=i 2^i\ e i\ 2 inducing a kernel U (k) as in (A5) of the Ap- 
pendix. 

There is another symmetry-allowed transitio n 13 ' 15 , namely 
the square/centered rectangle. In addition to strain, as in the 
square/rectangle, it also involves a shuffle because of the cen- 
ter site, so is not considered here. 



C. Rectangle/oblique case: Ny = 2, Nop = 1, n = 2 

The shear physical distortion e 3 changes a rectangle to an 
oblique polygon, and is the single Nop = 1 order parame- 
ter. The point group G =p2mm with vq = 2 elements goes 
to the subgroup g= p2 with v g = 1 elements, so there are 
Ny = 2 variants. The Jl{^) Landau part is the same as the 
square/rhombus case of (3.8); howeverthe n = 2 non-OP con- 
tributions are harmonic in the combinations e± = (e\±e%)/2, 
as fnon = ~2~ e + 2 + ~ir e - 2 - This yields a different compati- 
bility kernel^in f compat = ^U(k)\e 3 {k)\ 2 , as in (A6). 



IV. TRANSITIONS WITH N v = 3 

We consider two transitions with Ny = 3 variants, and 
Nop = 2 order parameter (OP) components, but with dif- 
ferent numbers n of non-OP strains. They are the (a) cu- 
bic/tetragonal transition in 3D (n = 4); and (b) trian- 
gle/centered rectangle in 2D (n = 1). See Figs 1 and 2. 



A. 3D cubic/tetragonal case: 

Nv = 3, Nop = 2,n = 4 

There are three axes along which the cubic unit cell can 
elongate, to make a tetragonal cell, so one expects three vari- 
ants, as in Fig 2. The cubic symmetry group G = Pmim with 
vq = 48 elements goes to the tetragonal group g = PA/mmm 
with v g = 16 elements, so there are& Ny = 3 variants, 
as expected. The variants are generated by joint action of 



the Nop = 2 order parameters that are the two 3D devi- 
atoric strain a 13 ' 18,28 , with the vector in OP space chosen as 
e = (ea, e 2 ) ~ (^{X 2 + Y 2 2Z 2 }, ^{X 2 - Y 2 }). 

The cubic/tetragonal Landau free energy has been consid- 
ered by Barsch and Krumhansl and other o 18 ' 28 . The invari- 
ants {I p } under the cubic point group, up to a maximum or- 
der p max = 4, are I 2 = e 2 + e\ = e 2 , I 4 = Jf , and 
a third-order invariant 73 = eg — 3e 2 e 3 . This is explic- 
itly seen to be a scalar under cubic-symmetries, as it can 
be written with (2.6a) in terms of invariants (XY Z) 2 and 
X q + Y q + Z q with g = 2,4, 6. From BK scalings as be- 
fore, and sign choices a 3 = —1,(74 = +1 f° r three minima, 
we have /l = (r — l)e* + /o(e3, e 2 ), where the temperature- 
independent fo of (2. 18a) is 

fo = h ~ C3I3 + C4/4, (4-1) 

and the unsealed elastic constants C^ 3 \ are related to the 
scaled ones as 

C 3 = C^\ 3 /E ; C 4 = C^\ 4 /E . (4.2) 

In polar coordinates in OP spaced, with e = ( 63,62) = 

e(cos 4>, sin (f>) of magnitude 

£^|e1 = [e 3 2 + e2 2 ] 1/2 , (4.3) 

the symmetry in OP space is manifestly carried by the third- 
order invariant, 

I 3 = e 3 (cos 3 4> - 3 cos 4> sin 2 0) = e 3 cos 3</>. (4.4) 
Then fo in the form of (2.23), with 773 = cos 30 is 

fo = e 2 - C 3 e 3 + C 4 e 4 + A/ ; (4.5a) 

A/o = C 3 e 3 (l-r /3 ). (4.56) 

The angular dependence is in A/o ~ — cos 3(f). The radial 
minima and degeneracy conditions on fo with rj 3 ((j> m ) = 1, 
yield 2 - 3C 3 + 4C 4 = and 1 - C 3 + C 4 = 0, fixing the 
two coefficients as 



C 3 = 2; 0-4 = 1, (4.6) 
achieved by choosing scaling parameters 

A = C^/2C^; E = C^(C {3) /2C^) 4 . (4.7) 

Then fL(e 3 , e 2 ) and hence e(r) is independent of material 
constants, i.e. there is quasi-universality of the first kind. The 
Ny = 3 variant minima of fo(e, cf>) are at angles 773 (<f> m ) = 
cos 3<f> m = 1> and radius e m = e(r), where 
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sin30 m =0; <j> m = v AT ' , m = 1, 2, 3(= N v ); (4.8) 



e(r) = |[l + Vl-8r/9].(4.9) 

For r = 1 , e(r) = 1, as required, and the upper spinodal is 
universal, t = t„ p = 9/8. The saddle-point barriers are at 
radius £b(r) = |[1 — ^/l — 8r/9], and angles 0f, Tn = (2m — 
1)tt/3. 

At transition the minima of the Landau free energy fall 
on the Ny + 1 = 4 vertices and center of a rightward- 
pointing equilateral triangl o 13,28 inscribed in a unit circle in 
Nop = 2 dimensional OP space, with corners at (e3, e 2 ) = 
(l,0),(-l/2,±\/3/2). (For a different sign choice ct 3 = +1, 
the triangle merely changes direction.) See Fig 3 and Ta- 
ble I. The upper bound on the possible types of domain wall 
between pairs of variants is Nw = Ny\/[2\(Ny — 1)!] = 
AV(AV-l)/2 = 3. 

The scaled Landau free energy ft = {r — l)e 2 + fo in OP 
components with the choice (4.6) is 

fa =I 2 - 2I 3 + I 4 = e 2 + e 2 - 2(e\ - 3e 3 e 2 ) + (e| + e 2 ) 2 

= (1 + 2e 3 )(3e 2 - e 2 + 2e 3 - 1) + (e 2 - l) 2 , (4.10a) 

where the second equation 13 explicitly shows the fourfold de- 
generate roots of /l = at transition. 
In polar coordinates as in (2.23), 

f L (s, 0) = [(t - l)e 2 + e 2 (e - l) 2 ] + A/ 0) (4.105) 

where A/ = f (e, (f) - / (e, </> m ) = 2(1 - r; 3 )e 3 vanishes in 
the angular directions = cf> m of the minima. As fz, (e, <f>) is 
material-independent, there is again quasi-universality of the 
first kind. Here from the definition of (2.22d), gp = r — 1 + 
(e-1) 2 . 

The Ginzburg term is f G = £ 2 [(Ae 3 ) 2 + (Ae 2 ) 2 ]. There 
are n = 4 non-OP compressional e\, and shear strains 
64, es, ee, that can be written in terms of the OP strains, so 
e « = Tlt=2 3 Ba e t w i m ^ = 1,4:, 5, 6. Substituting into the 
non-OP harmonic terms yields the cubic/ tetragonal potential 
with the 2 x 2 matrix kernel of (2.19), given in (A23) of the 
Appendix. The 3D relaxational OP strain simulations-^ did 
not explicitly state the kernel, now given here for complete- 
ness. 



shear strains, with the OP vector chosen as e = {e%, e 3 ) ~ 
\{X 2 — Y 2 ,XY). The triangle point group G= p6mm with 
vq = 6 elements goes to the centred rectangle subgroup g = 
c2mm with v g = 2 elements, so there arei^ Ny = 3 vari- 
ants, as expected. The third-order invariant under ir/3 rota- 
tions of the triangular lattice is now J 3 = e\ — 3e2e|. For 
(X,Y) = R(cosa, sin a), one finds J 3 = (i? 3 /8) cos 6a, 
manifestly invariant under the a —> a + ir/3 triangular sym- 
metry. The scaling parameter choices are as before. The fi- 
nal scaled Landau free energy^ 3 - is formally similar to the cu- 
bic/tetragonal case of (4.10a), with e 2 and e 3 interchanged, 

fo = h~2h+h = el+el-2{el-Ze 2 el)+{e 2 2 +el) 2 . (4.11) 

At transition the Nv + 1 = 4 degenerate Landau minima 
again fall on the three vertices and at the center of an equilat- 
eral triangle inscribed in a unit circle. The maximum number 
of domain wall types is Nw = 3. 

The Ginzburg terms are the same as above, while the 
compatibility potential in 2D, is of course different. The 
single n = 1 non-OP (compressional) strain is e\ = 
— J2e=2 3 OieijO\ from the 2D compatibility constraint. 
Substitution into the harmonic term as in (2.19) immediately 
yields^, f non = (Ai/2)\ei(k)\ 2 -> f com pat(e 2 ,e 3 ) with a 
2x2 compatibility-kernel matrix of (A2). 

V. TRANSITIONS WITH N v = 4 

We consider two transitions with Nv = 4 variants, and dif- 
ferent order parameter (OP) components. They are the (a) 
cubic/trigonal case in 3D ( Nop = 3, n = 3); and (b) 
square/oblique case in 2D (Nop = 2, n = 1). See Figs 1 
and 2. 



A. Cubic/trigonal case: Nv = 4, Nop — 3, n = 3 

The distortion acts along body diagonals of the cube and the 
Nop = 3 shears are the three components of the OP vector 
e =_(e 4 ,e 5 ,e 6 ) - (YZ,ZX,XY). With the cubic group G = 
Pm3m with v G = 48 going to the trigonal or rhombohedral 
subgroup g = P31m with v g = 12, there ara& Ny = 4 shear- 
induced variants. The four invariants up to order p max = 4 
are I 2 = e\ + e| + e§; h = e4e 5 e 6 ; I 4 = e| + ef + e\\ and 
l' A = I 2 2 , with N mat = 4 material coefficients. The scaled 
Landau free energy is 

h = (r- 1)I 2 + fo(e 4 , e 5 ,e 6 ), (5.1a) 

where 



B. Triangle/centred-rectangle case: Nv = 3, Nop = 2, n = 1 

There are three ways to convert an equilateral to an isosce- 
les triangle, with the unit cell of the (equilateral) triangle 
changing to a centered rectangle, so one expects three vari- 
ants, as in Fig 1. The OP are 2D Nop — 2 deviatoric and 



fo =h~ C 3 I 3 + C' 4 ll + C4I4, (5.16) 
and the scaled parameters are related to unsealed ones by 

C(3) „ f7'(4) (7(4) 

C 3 = A -— - ; C 4 = A — = — ; C 4 = X ——. (5.2) 

-CyQ £jQ -C'O 
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In spherical polar coordinates in OP space, 

e= (e4,e5,ee) = e(sin9 cos <j>, sin sin 0, cos 6), (5.3a) 
with magnitude 

e=\e\ = [e 4 2 +e 5 2 + e 6 2 } 1/2 . (5.36) 
The invariants in OP space are then 

I2 = £ 2 ; I3 = e 3 sin 2 9 cos 9 sin <f> cos 0; I' 4 = e 4 , 

h = e 4 [sin 4 6»(cos 4 + sin 4 0) + cos 4 6\. (5.4) 
Using trigonometric identities, /o is 



C 

fo(e,9,4>) =£ 2 + C' 4 e i - -^e 3 sin (9 sin 26> sin 20 



+C 4 e 4 [l - i sin 2 20 - - sin 4 9 sin 2 20], (5.5) 

and has three remaining material constants C3, C4, C 4 , while 
there are two scaling parameters A,Eq. 

The radial minimum dfo/de + 2(r — l)e = 0, and r = 1 
degeneracy condition /o(e = 1, 9 m , rn ) = yield 



C 3 = 6^3; C 4 ' = 1 - C4/3, (5.6) 
achieved by the choice of scaling parameters 



A 



jf ^ . Eo . A 3^. (5 . 7) . 



The angular minimum conditions dfo/d(f> = 0, and 
dfo /d6 = yield minima at cos 20 m = or at 



have equal relative separations or polyhedral sides of e.g. 
I el — I = 2^/2/3 e m , and relative cosines of eg el.ei/e^ = 
cos?/; = —1/3, so the angle between vectors is the well- 
known tetrahedral angle ip = eos _1 (-l/3) = 109°28'. 

At transition, the Landau minima fall on the vertices and 
center of a tetrahedron, inscribed in a sphere of unit radius 
in Nop = 3 dimensional order parameter space, as in the 
schematic points of Fig 3. The generic numbers for transition 
are in Table I. The maximum number of domain wall types 
is Nw = 6. Of course in this and other cases, this is only an 
upper bound, and energetic considerations might lead to fewer 
types actually appearing in the final microstructure. 

Then the scaled Landau free energy in components is 

fh = (t-1)I 2 +{I2-6V3I 3 + (1-C 4 /3)I 2 +C 4 I 4 }. (5.10a) 
In polar coordinates as in (2.23), it is 

h = [(t - l)e 2 + e 2 (e - l) 2 ] + A/ , (5.106) 

manifestly showing the degeneracy at transition 
/i(l,0m,^m) = 0. The angular part A/ = 

fo(e,4>,0) - fo(e,<f) m ,6 m ), is 



A/ 



3^3 3r 



3V3 



+C 4 s 4 



sin 2 26> 



sin 9 sin 29 sin 20] 



■sin 4 9 sin 2 20]. (5.11) 



This carries material dependence through C4, but vanishes at 
minima, so /l(£, 9 m: <f) m ) and hence e(r) are still universal 
for all r i.e. there is quasi-universality of the second kind. 
Here from (2.22d), g L = r - 1 + (e - l) 2 . 

The material-independent OP magnitude at the tetrahedron 
corners is 



sin' 



= 0, 3cos 2 m = l, (5.8a) 



4> m = (2m - l)ir/N v , m = 1,2,3,4(= N v ). (5.86) 

For appropriate positive signs of the second derivatives, there 
are two points each on the northern and southern unit hemi- 
spheres in OP space, 



e(r) = -[l + Vl-8r/9], (5.12) 

and happens to be the same form as for the cu- 
bic/tetragonal transition, with the saddle-point barrier at 
£~b(r) = f [1 - y/1 - 8r/9]. The Ginzburg term is f G = 
£ 2 [J2i=4 5 e{^- e i) 2 ]- The n = 3 non-OP strains are the com- 



pressional (ei) and the deviatoric strains (e2,e3). The 3x3 
compatibility kernel Uw(k) of (2.19) is given in (A20). 



m = tt/4, 5tt/4, 0„ 



B. Square/oblique case : N v = 4, Nop — 2, n — 1 



4> m = 3tt/4, 77t/4, 9 m =9 + tt, (5.9) 

where cos 9 = l/v3. These four minima m = 1,2,3,4 
are at the ends of vectors in three-dimensional OP space 

e = e m = £ m (sm0 m cos0 m ,sin0 m sin0 m ,cos0 m ), that 



A square converts to an unequal-sided oblique polygon with 
the deviatoric and shear e 3 order parameters acting simul- 
taneously. Since each can distort the square in two ways, one 
expects four variants, as in Fig 1 . The harmonic parts are as- 
sumed to soften at the same T c (otherwise they would be two 
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separate, unrelated transitions). The joint action of the two 
order parameters is through their coupling. The square point 
group G= pAmm with vq = 4 goes to g = p2 with v g = 1, 
so there are 1 ^ Ny = 4 variants, as expected. The anharmonic 
invariants are as in the square/rectangle case, e% , e 2 and ef, 
e| or pmax = 6 separately for each OP. We consider the sim- 
plest case of equal elastic constants, and the simplest coupling 
C2 2 e§. The scaled free energy term fo is: 

fo = (el+el)-C 4 (ei+el)+C\(e 6 2 +el)}-C' 4 e 2 2 e 3 2 , (5.13) 

where the scaled and unsealed coefficients are related by C4 = 

CMA 4 /£ ; C 6 = C^X 6 /E ; C 4 = C' {i) X 4 / E Q . 

Transforming to polar coordinates (e2,e 3 ) = 
e(cos </>, sin </>), with r/ 2 = cos 20, and using trigonometric 
identities, this becomes 

fo =t a {s)+U(e) m 2 , (5.14) 
where the coefficients are 

t a = s 2 -^(2C i + C' i ) + ^C 6 , (5.15a) 



U = ^C 6 e . (5.156) 

The angular dependence is /o ~ cos 4</>. The degeneracy 
condition, and the radial minimum condition at transition fi- 
nally yield C 4 + CJ2 = 4; C 6 = 4, achieved through the 
choice 

A 2 = (C (4) +C' (4) /2)/C (6) ; E = A 6 C (6) /4. (5.16) 
The angular minima are at i]2{4>m) 2 = 1 or 



sin# m = 0; <f> m = ± — — -L-, m = 1, 2, 3,4(= Ny). (5.17) 
i\y 

Then the scaled free energy in OP components is 

fo = S 2 -(4-^C' 4 )(el+et)+i(e 6 2 +e 6 3 )~C' 4 e 2 2 e 3 2 . (5.18a) 
In polar coordinates as in (2.23), 

/ L = [(r-l)e 2 + £ 2 ( £ 2 -l) 2 ] + A/ , 

Af = e 4 (3s 2 -2 + C'J2) cos 2 2cf>, (5.186) 

so the quasi-universality is of the second kind, with /z,(e, 4> m ) 
and e(r) of (3.6), both independent of material constants. 
Here g L = r - 1 + (e 2 - l) 2 . 



At transition, the Ny = 4 variant minima fall on the ver- 
tices of a square inscribed in a unit circle in Nqp = 2 di- 
mensional order parameter space. These are the only vari- 
ants for C4 > 1, as the trivial roots on the axes (e2, 63) = 
(±£, 0), (0, ±e) are then unstable. The maximum possible 
number of domain wall types is Nw = Ny(N\/ — l)/2 = 6. 

The Ginzburg term is fo — £ 2 [(Ae2) 2 + (Aea) 2 ]. There is 
a single n = 1 non-OP variable ei as in the triangle-centered 
rectangle case, so the compatibility kernel from the compres- 
sional harmonic term e\ is the same as the 2x2 matrix Uw 
of (A2). 



VI. TRANSITIONS WITH N v = 6 

We consider two transitions with Ny = 6 variants, and 
No p = 2 order parameters (OP), but with different numbers 
n of non-OP strains. They are the (a) cubic/orthorhombic case 
in 3D (n = 4); and (b) triangle/center rectangle case in 2D 
(n = 1). See Figs 1 and 2. 



A. Cubic/orthorhombic case Nv = 6, Nop — 2, n = 4 

For a cubic to orthorhombic distortion, the cross-sectional 
area perpendicular to each of three axes can be rectangular 
in two ways, so one expects six variants, as in Fig 2. The 
symmetry group G = Pm3m with vg = 48 elements goes to 
the orthorhombic group g = Pmmm with v g = 8 elements, 
so£ Ny — 48/8 = 6, as expected. The variants are gener- 
ated by combined action of the Nop = 2 order parameter 
components that are the two 3D deviatoric strains, with the 
vector in OP space chosen as in the cubic/ tetragonal case, 
e= (e 3) e 2 ) ~ (^{X 2 + Y 2 - 2Z 2 }, ^{X 2 - Y 2 }). 

The cubic/orthorhombic free energy in Cartesian strains 
has been considered for fitting to FePd experiments^; here 
however, we work with physical strains. The previous cu- 
bic/tetragonal case of (4.1), with a third order invariant ^3 = 
e 3 a — 3e 3 e2 2 , yielded three minima alternating with three 
maxima, on the unit circle. For six minima on the unit cir- 
cle, a sixth order invariant 1$ = if will be the leading angular 
term. 

We consider two cases, with up to sixth order, and up to 
eighth order invariants. For invariants of up to p max = 6th 
order, the free energy has N max = 4 material coefficients. 

Fl = (t- l)(e§ + e|) + / (e 3 , e 2 ) where f is 

fo = h- CJZ - C 6 I 6 + C' 6 ll (6.1a) 

with signs (74 = ere = — 1, a' 6 = +1. For materials with other 
coefficient signs, we are forced to go to higher p max = 8th or- 
der, and the additional invariants are 1$ = I 2 I 2 , 7g = 7| ■ (The 
odd invariants 13,15 = I 3 I-2,h = I5I2 give sign-varying 
contributions to derivatives dfo / d<fi at different minima that 
should be equivalent, so we set their coefficients to zero from 
the start.) 
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The scaled free energy up to eighth order, with N max = 6 
material coefficients is 

fo = h- Cdl + C 6 I 6 - C' 6 li - C 8 I 8 + C' & lt (6.1b) 

It is convenient to define Co = C 6 — Cq,Cs = 
Cg — Cg. Transforming to polar coordinates, e = (e3, e^) = 
e(cos 0, sin 0) we get I 3 = e 2 cos 30 as before, so 7 6 = 
£ 6 ?73, = £ 8 f?l where 773 = cos 30. 

Collecting terms, the sixth and eighth order cases of (6.1a) 
and (6.1b) can both be written as 

fo = t (e)+t 6 (e)(l- V 2 3 ), (6.2) 

where from (6.1a) 

to = e 2 - C 4 e 4 + C 6 (_) £ 6 ; t 6 = C 6 e 6 ; (6.3a) 

while from (6.1b), 

t = s 2 -C 4 e 4 -Ce { ~ ) s 6 +C 8 i ~ ) e 8 ; t 6 = C s e 8 -C 6 e e . (6.36) 

In both cases, the angular dependence is fo ~ — cos 60. 

For the sixth order case of (6.3a), the degeneracy fo(s = 
1,0 = 0m) = and radial minimum condition dfo/de + 
2(t — l)e = determine two of the constants as Cq^ = 
1, C4 = 2. The scaling parameters for these values are, with 
unsealed elastic-coefficient ratio a = C (4) /(C" (6) - C (6) ), 

A 2 = q-/2, E Q = X 4 C {4) /2. (6.4) 

Then the angular contributions clearly yield dfo /<90 = 
roots at 

sin60 m = 0; m = 2 ( m m = 1) 6 (-_ (6.5) 

ivy 

where 77 3 (0m) 2 = 1. 

The scaled free energy in components in OP space is, with 
the above scaled coefficients, 

f L = (r - l)I 2 + I 2 {h - l) 2 + C 6 (/ 2 3 - J 3 2 ). (6.6) 
In polar coordinates as in (2.23), 

/ L = [(T-l)e 2 + £ V-l) 2 ] + A/ , 



A/o = -C 6 e 6 (l-cos60), (6.7) 

The last term A/o vanishes at the six minimal directions, 
where the material constant Co is eliminated, so there is quasi- 
universality of the second kind. The OP magnitude e(r) is 
as in the tetragonal/orthorhombic case of (3.6). The variants 



are manifestly degenerate with austenite as /l(1, 4>m) = at 
transition. 

For the eighth order case of (6.3b), there are three mate- 
rial constants C4, C^~\ C^~' in /z,(e, 0m )> an d only two 
remaining scaling parameters A, Eq. The degeneracy and ra- 
dial minimum condition now determine two of the constants 
as Co { - ] = 3 - 2C 4 ,C* 8 ( ~ ) = 2 - C 4 . The equivalent 
condition 3C 8 (_) - 2C* 6 ( ~ ) - C 4 = 0, yields with (2.12), 
a quadratic, 3A 4 — 27A 2 — ory, where a is as above and 

7 = (C' (6) - C^)/{C' [S) - C^). The positivity of sec- 
ond derivatives requires that C% > Co and 2 — C4 > while 
A 2 > below, further requires 3 — 2C4 > 0. 

However A(<C 1) can also be obtained from the relation 
between scaled and unsealed coefficients (2.12), as A 2 = 
a(3 - 2C 4 )/C* 4 = 7(2 - C* 4 )/(3 - 2C 4 ). Demanding consis- 
tency yields A, Eq, C 4 in terms of the unsealed elastic coeffi- 
cients, but here C 4 is no longer just a universal number. The 
scaling parameters are then 

A=( 7 /3)[l + {l + 3a/ 7 } 1 / 2 ]; E = C^X 4 /C 4 ; (6.8a) 

C 4 = (3/2)/[l + (7/6o)(l + {1 + 3a/ 7 } 1/2 )]- (6.86) 

As elastic constants vary, the constant C 4 (a/7) moves in a 
narrow range 3/2 > C 4 > 0, e.g. C 4 = 1 for a/j = 1. 

The scaled free energy in OP components for the eighth 
order case from (6.3b) is then 

h = (j- l)h + h{h - 1) 2 {1 + (2 - Ci)I 2 } 

+ (C S l2-Co)(I 2 3 -h 2 )- (6.9) 
In polar coordinates as in (2.23), 

h = [(t - l)e 2 + e 2 (e 2 - 1) 2 {1 + (2 - C 4 )e 2 }} + Af , 



A/o = ^(C 8 e s - C 6 e 6 )(l - cos 60). (6.10) 

with Cs, Co eliminated at minimal angular directions. 
The root e(r) is the solution of a cubic 

4(2-C 4 )X 3 +3(5-2C 4 )X 2 +2(3-C* 4 )A:+(t-1) = 0, (6.11) 

where X = s 2 — 1 > 0. At r = 1, the required root is X = 0, 
and just below transition is linear, X ~ (1 — t)/2(3 — C 4 ). 
Close to zero temperature, with t = — |r| and for |r(T = 
0)| = T C /(T - T c ) » 1, one has X ~ [|r|/4(2 - C4)] 1 / 3 . 
Thus the scaling procedure carries through even for the eighth 
order case, but there is now quasi-universality of the third 
kind, with a weak residual material-dependence through C 4 . 

At transition and for both cases, the Ny = 6 nontrivial 
Landau minima fall on the vertices of a hexagon^ inscribed in 
a unit sphere in Nop = 2 order-parameter space as in Fig. 
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3 and Table I. The upper bound on the number of possible 
martensite-martensite domain wall types is Nw = Ny(Ny — 
l)/2 = 15, although as mentioned, not all of these may be 
seen, for energetic reasons. 

The Ginzburg term, and the compatibility potential from 
the n = 4 non-OP harmonic term is the same as in the cu- 
bic/tetragonal case of (A23) of the Appendix. 

B. Triangle/oblique case Ny = 6, Nop = 2, n = 1 

The two order parameter components are the single devia- 
toric and single shear strains, e = (e2,es) = e(cos <fi, sin </>). 
The triangle point group G = p6mm with vq = 6 goes to the 
subgroup g =p2 with v g = 2 so there are-^ Ny = vojv g = 6 
variants. (The ratio of the numbers of rotational elements 
would however give 3/1 = 3 variants.) 

The scaled Landau free energy density^ 3 , is formally simi- 
lar to the cubic/orthorhombic case of (6.7). At transition the 
Ny + 1 = 6 + 1 degenerate Landau minima again fall on 
the six vertices and at the center of a hexagon inscribed in a 
unit circle, as in Fig 2. The maximum number of domain wall 
types is N\y = 15, but not all may finally appear. 

The Ginzburg term, and the compatibility potential from 
the n = 1 non-OP harmonic term is the same as in the trian- 
gle/centered rectangle case of (A2). 

VII. PSEUDOSPIN HAMILTONIANS 

The idea of using discrete-variable pseudospins to approx- 
imate continuous-variable distortions on a lattice was pro- 
posed earlie r 20,21 , and has been pursued 2 ^—. We had sug- 
gested obtaining pseudospin hamiltonians for various transi- 
tions by substituting the order parameter (OP) values at the 
polar coordinate Landau minima, into the total scaled free 
energy 22 . The Nop = 1 case for the square/rectangle case 
as e2(f) — > sS(r), where the three minima induce a spin-1 
model S = 0, +1, —1 with three values on a line, of a single- 
component pseudospin. With S 6 = S 4 = S 2 = 1 or 0, 
the nonlinearities in the Landau free energy collapse to — > 
z( t ) 2 9l{t)S 2 where <7l(t) = t — 1 + (e 2 — l) 2 changes sign 
at transition. The Ginzburg and compatibility terms are also 
written in terms of pseudospins. This yields a temperature- 
dependent pseudospin hamiltonian i/( 5) = F(e.2 — > SS) that 
is like a generalized Blume-Capel spin-1 model 19 . The hamil- 
tonian has a temperature-dependent quadratic term, a nearest- 
neighbor ferromagnetic term, and a PLA compatibility term 23 , 

r 

+ J2^-U(r-nS(r)S(n], (7-1) 

r,r' 

with gradient V realized as difference operators A on 
a computational grid, as mentioned. Here Dq(T) = 



2e 2 (T)Eo/kBT, and the hamiltonian is diagonal in Fourier 
space, 

l3H=^Y,Qo(k)\S(k)\ 2 ; (7.2a) 

k 



Q (k) = D [g L (r) + ?K 2 + -j-U(k)]. (7.2b) 

Here on a grid of unit lattice constant, Ku = 2sin(fc p /2), 
with n = x,y. This hamiltonian has been studied in a lo- 
cal meanfield approximation, under a cooling ramp obtaining 
glassy domain-wall textures, dependent on cooling rate and 
initial conditions^ 3 -. 

For No p > 1 , Ny > 2 transitions we do not simply get a 
generalized spin-j model with 2j + 1 states on a line, where 
j = Ny/2. Instead we obtain clock-like models^ with dis- 
crete S vector variables pointing to Ny + 1 corners and centre 
of a polyhedron in A^p-dimensional space, as denoted by ar- 
rows in Fig. 3. Since the zero state is included, these may be 
termed 'clock-zero' Zjv v +i models. Note that, unlike pure 
clock Zjv models, the spin square-magnitude S 2 (r) is still 
a statistical variable and not a constant, because of the zero 
states. Choosing Nop = 2, 3 component strains only at the 
minima induces vector pseudospins in OP space, 

e(r)^eS(r). (7.3) 

The variant angular dependence fr, ~ — cos Ny(f> generates 
the clock- variable directions. 

The general temperature-dependent pseudospin hamilto- 
nian is 



H(S l (r))=F(e t ^eSi). (7.4) 

As in the square/rectangle case, the radial part of the Landau 
term with OP nonlinearities, collapses to a quadratic in the 
pseudospin magnitude, since nth powers of the spin-vector 
magnitude \S\ n = \S\ 2 = 0, 1. 

Although in zero stress the uniform state is no longer a Lan- 
dau minimum below the lower spinodal T c , there is a possibil- 
ity that nonuniform textures exert local internal stresses to fa- 
vor the zero value at a site even at low temperatures. Also, the 
original free energy in OP strain always has a turning point at 
the origin to support dynamical transient zeros, that although 
few in number, could play a catalytic role in microstructural 
evolutions 3 -. Hence we retain zero spin values at all temper- 
atures, allowing their permanent/transient existence to be de- 
termined dynamically. 

The hamiltonian in coordinate space is 

PH=^[Y J {9L(r)S e (r) 2 + e(AS e ) 2 } 
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',f,i,i' 



(7.5a) 



and is transition-specific through the S values, the temper- 
ature dependence of <7l(t), and the compatibility potential. 

Note that the anisotropic terms S(r).u (f—r').S(r l ) are from 
compatibility anisotropies in OP space, and differ from mod- 
els with electric dipoles d that have anisotropies relative to co- 
ordinate space axes 2 ^ {d(r).f}{d(r / ).r'} /\r — r'\ d , (although 
the powerlaw fall-offs with exponent d, are the same). 
The hamiltonian is diagonal in Fourier space, 



h,u>(k)St(k)S e ,(ky; (7.6a) 



At. 



QoMQ = A)[{fli(T)+W^ + -fUw(k)]. (7.66) 

Transitions with free energy quasi-universality of the first and 
second kind, have reduced pseudospin hamiltonians with uni- 
versal coefficients g^ of the on-site term. Apart from the over- 
all Eq, the material-dependence is only through the texture- 
inducing Ginzburg and St Venant terms. 

For two-component OP cases in addition to S = 0, the 
pseudospin is S = (cos cj) m , sin (j> m ), where m takes on Ny 
values. Thus for the three variants of the cubic/tetragonal 
and triangle/center-rectangle transition, there are three val- 
ues cj} m = 0, 2ir/3, 47r/3 on corners of a triangle, and ol = 
r— l+(e— l) 2 where e(r) is from (4.9). The four spin vectors 

are£=(0,0),(l,0),(-±,±^). 

For the square/oblique transition there are four values 
4> m = 7r/4, 37r/4, 57r/4, 77r/4 on corners of a square, and 
cjl = t — 1 + (e 2 — l) 2 , with e(r) of (3.6). The five spin 

vectors are (0,0),(±^,^),(±^,-^). 

For the three-component OP of the cubic/trigonal tran- 
sition, the pseudospin vectors are non-planar as S = 
(sm9 m cos4> m ,sm0 m sm<f) m ,cos9 m ), where 9 m ,<j) m take 
on four values of (5.8) at the corners of a tetrahedron. The 
coefficient = r — 1 + (e — l) 2 , with e(r) the same as 
the cubic/tetragonal case of (4.9). The five spin vectors are 



(0,0,0),(± : 



I 111 ,/5> T . 



V5'" L v^ ! V3 ,,l " L ^' T v'3 r " V^' 

For the two-component OP and six variants of the cu- 
bic/orthorhombic and triangle/oblique transitions, <f> m = 
0, 7r/6, 27r/6, 57r/6on the six corners of a hexagon, and the 

seven spin vectors are (0, 0), (±1, 0), (±, ±^), (-§, ±^), 
where for the p max = 6 case, <7l(t) = r — 1 + (e 2 — l) 2 is 
universal, and e(r) as in (3.6). For the p m ax = 8 case, with 
quasi-universality of the third kind, there is a weak material- 
dependence through g L = r - 1 + (e 2 - 1) 2 [1 + e 2 (2 - C 4 )] 
with 3/2 > C 4 > 0, and e(r) as in (6.11). 

The pseudospin reduced hamiltonians here are not just writ- 
ten down, but are induced from the scaled free energy, that 
encodes the specific symmetry, nonlinearity, and compatibil- 
ity of each ferroelastic transition: a continuum-variable mate- 
rials science model is mapped to a discrete-variable statistical 
mechanics hamiltonian. 



VIII. SIMULATIONS OF STRAIN TEXTURES 

Spatially varying strain textures in ferroelastics can be nu- 
merically simulated in continuous strains by free energy re- 
laxations, or by discrete-strain pseudospin hamiltonians. 



A. Scaled free energy relaxations 

Of course, there has been much simulation work in the dis- 
placement representation or with phase fields 2 ^, and one com- 
mon problem is the choice of the many N mat material coef- 
ficients, that require fitting to experiment for each material 7 . 
With scaled free energy strain dynamics or with pseudospin 
hamiltonians, the material-dependence is essentially elimi- 
nated or reduced, to group many materials with the same 
structural transition in the same quasi-universality class. Thus 
one does not have to explore the full many-parameter space 
of N mat unsealed coefficients, and simulations by different 
groups can be more easily compared. 

Local equilibrium microstructures can be found from 
relaxational dynamics of unsealed order parameter (OP) 
strains, {e^(r, t)} and the unsealed free energies: de' e /dt — 
-T^dF/de' e , where is a kinetic constant of dimension 
inverse-energy x time. Scaling strains, e' e = \ei and energies 
F = EqF, where F = Fl+ Fq + F compat , we obtain a char- 
acteristic decay rate To = r'°) Eq / A 2 . Then with a dimen- 
sionless time t = Tot absorbing Eq, a scaled dimensionless 
dynamics is obtained, 



de e 
dt 



dF{e e ) 
dee 



The underdamped dynamics- 3 could be similarly scaled to 
quasi-universal, dimensionless form. 



B. Pseudospin simulations 

Microstructures can be studied using pseudospin hamilto- 
nians, by solving self-consistency equations from local mean- 
field approximations and by Monte Carlo spin simulations^. 



1. Local meanfield 

The completely uniform k = meanfield contribution to 
the hamiltonian is H ~ f L (t)S 2 (k = 0), as the Ginzburg and 
compatibility terms vanish. Thus as Jl = £ 2 gL changes sign 
at transition, we have < S 2 >= above the transition, 

and < S 2 >= 1 below the transition, faithfully reproduc- 
ing the displacive transition of the original strain variable 23 . 
Strain textures with k ^ 0, can be captured by local mean- 
field cr^(r) = < Si(r) > approximations in both coordinate 
and Fourier space. Each spin sees a local meanfield, and there 
is a subtraction for consistency of averages < SS >~ era: 
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-az(r)<j e ,(7) (8.2a) 
S e (k)Se(k)* -> S t (k)ai>{k)* + ai{k)St>{k)* 

-ai{k)a e {ky. (8.2b) 

This is equivalent to substituting 5 = <r + (55 into the hamil- 
tonian, and linearizing in SS = S — a. 
The meanfield hamiltonian is then, 

PH MF = Qt(f)S t (f) = Y, Qt(k)*S t (k), (8.3) 

if e,k 

where an additive constant — | Qi(r)o-g(r) on the right is 
suppressed. Here Q t (r) = Y.e,^ Qo,w(r ~ f^t'ir 1 ), or 
Qi(k) = J2e> Qo,i£'{k)(Je>{k). The same results are obtained 
from a k ^ variational approach in the weights p(k), mini- 
mizing Trace(H p) — J2k PW ^ n P(k)> mat nas avera g es sucn 
as cr^(fc) = Trace{p(k)S(k)}. 

The self-consistency condition for the local, nonuniform 
meanfield is then 

For the AVjip = 1, jVy + 1 = 3 square/rectangle case, one 
obtains^ 

cr(f) = -2sinhQ(f)/[l + 2 cosh Q(f)]. (8.5). 

2. Monte Carlo simulations 

Metropolis algorithms for nearest-neighbor Ising models 
with spin components No p = 1 and values Ny = 2, are pro- 
totypical. The Monte Carlo method can be applied to spins of 
Nop = 1,2,3 components, with Ny + 1 = 3, 4, 5, 7 values 
and powerlaw anisotropic interactions. Fast Fourier Trans- 
forms (FFT) are used to easily treat the full PLA compatibil- 
ity potentials; this is preferable to attempting an uncontrolled 
truncation of the 1/R d interactions to some arbitrarily far- 
neighbor interaction, while staying only in coordinate space. 
The procedure is: 

(i) Flips of an input spin configuration with an input energy 
Hi n p U t are made on the lattice in coordinate space, with equal 
probability over all Ny + 1 values at a site, yielding a trial 
new configuration {S(r}}. 

(ii) An FFT gives the Fourier space {S(k)}. 



(iii) The trial energy H tr iai (diagonal in Fourier space as 
in (7.6b)) is then evaluated. The energy difference AH = 
Htriai — Hi npu t between configurations in the Boltzmann- 
factor e~P AH , determines the accepted configuration. The 
cycle then repeats. 

The sites are labelled with random numbers and arranged 
in an increasing-value chain, so sequentially visiting every site 
in the chain means every site is visited randomly, but once and 
only once, in 1 MC sweep, that then corresponds to 1 MC step 
per spin. 



IX. SUMMARY 

Generalizing a procedure due to Barsch and Krumhansl, we 
have shown that ferroelastic free energies can be scaled in di- 
mensionless form, defining a quasi-universality class for all 
materials with the same group/subgroup structural transition. 
Whereas unsealed free energies with fitted coefficients are 
specific to the fitted material^, the scaled dimensionless free 
energies are in a form relevant for all materials with the same 
ferroelastic transition. A simplifying approximation works to 
leading order in the (unsealed) spontaneous- strain magnitude 
at transition A, that is typically a few percent. To this order, the 
scaled Landau free energy minima turn out in most cases to be 
material-independent, depending only on the dimensionality 
d, number of order parameter components Nop', and num- 
ber of low-temperature structural variants Ny- The minima 
in order-parameter space fall on the corners of 'polyhedra' 
inscribed in 'spheres' in Nop dimensions, with radius the 
scaled order parameter magnitude, that is unity at transition. 
The scaled variational free energies in terms of the local order- 
parameter strain components, have Landau, Ginzburg, and St 
Venant powerlaw anisotropic interactions, and can be used in 
relaxational or underdamped dynamic simulations. The com- 
patibility kernels are calculated for all transitions considered, 
using a constraint-substitution method. 

The polyhedral arrangement of minima in the strain vari- 
ables immediately suggests a reduced description with strain 
free energies inducing 'clock-zero' models, with Nop- 
dimensional discrete pseudospin vectors at each site, point- 
ing to Ny + 1 possible states. The discrete-variable pseu- 
dospins in local meanfield and Monte Carlo simulations 
can be shown^i to reproduce previous textures obtained in 
continuous-variable dynamics- ^ 14 ] 28 . 

Further work could include pseudospin clock-zero simula- 
tions of more of the transitions considered here; the addition 
of quenched disorder; a determination 7 '^, of the scaled free 
energies and compatibility kernels for more of the 94 possible 
ferroelastic transitions 6 ; scaling of experimental results for 
different materials to explore data clustering; and simulating 
complex oxide models, through couplings of strains to other 
fields^. 

This work was supported in part by the U.S. Department of 
Energy, and by ICTP, Trieste. 
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X. APPENDIX: ST VENANT COMPATIBILITY KERNELS 
FOR FERROELASTIC TRANSITIONS 

We derive in this Appendix, the Fourier space compatibil- 
ity kernels, for increasing numbers n of non-order parameter 
strains: n = 1, 2 for 2D transitions, and n = 3, 4, 5 for 3D 
transitions. 

The background for compatibility ideas includes the fol- 
lowing. Kartha et at^ performed Monte Carlo simu- 
lations in the displacement-vector representation, using a 
square/rectangle variational free energy F(Vu) that, in ad- 
dition to a Landau term nonlinear in the order parameters 
(OP), also included harmonic terms in compression and shear 
~ (V x u x + VyUy) 2 and ~ (S7 y u x + V ' x u y ) 2 . They found 
diagonal domain walls. To understand the orientation, they 
followed Baus and LovetUS and went over to the strain rep- 
resentation, analytically minimizing these non-OP ~ ei 2 , e 3 2 
terms subject to the St Venant constraint using Lagrange mul- 
tipliers. The resultant square/ rectangle compatibility kernel 
explained the ir/4 orientation preference. However, strain- 
representation simulations in the OP e 2 , for a free energy 
F(e 2 ) that explicitly included the St Venant term, were not 
pursued. This changeover to the OP-strain working vari- 
able was done in relaxational simulations —, yielding the same 
diagonal domains. Strain-variable simulations using com- 
patibility kernels for other 2D case a 13 ' 15 , and the 3D cu- 
bic/tetragonal cas&ii, found textures as obtained in the dis- 
placement and phase-field representations 2 ^; and an under- 
damped strain dynamics including compatibility potentials 
was proposed, including Langevin dynamics noise terms with 
powerlaw spatial correlations^ 3 -. 

Here we pursue this strain-representation project, contribut- 
ing to a catalog of 3D ferroelastic No p x No p matrix kernels 
that incorporate the three compatibility constraints (2D ker- 
nels for ferroelastic transitions were previously given ini£). 
Instead of using, and solving for, three Lagrange multipliers, 
we directly solve the constraints in Fourier space for three 
non-OP strains, and substitute in the harmonic non-OP free 
energies, followed by an unconstrained minimization in any 
remaining non-OP strains 17 . The four 3D transitions have 
three distinct kernels, that are plotted in Fig 4. As a check, we 
apply the direct substitution method to find 2D kernels, previ- 
ously obtained by Lagrange multipliers^—. In Fourier space 
derivatives go as eg V M 2 — > — fc^ 2 , while on a grid with dif- 
ference operators A M 2 — > —K 2 where = 2sin(fc M /2). 
Here in compatibility equations we write for simplicity just 
the wave vectors like with the understanding that they can 
be replaced by in grid simulations. 

The compatibility constraint in minimizing the harmonic 

non-OP terms f non = X^i(^j/2)l e i(^)| 2 onr y affects nonuni- 
form strains with nonzero wave-vector. The uniform or zero 
wavevector parabolic terms are freely minimized by zero 
values, a(k = 0) = 0. Thus a prefactor of v{k) = 1 — St 
must be inserted in the results below for the kernels. The re- 
sultant sign-varying compatibility potentials have zero spatial 



average < U(R) >~ U(k = 0) - v(k = 0) = 0, 
rather than a divergence as for isotropic potentials 

<u(R) >~£^i/.R d ~iniV. 



A. 2D transitions 

In 2D, the Fourier constraint of (2.3b) in terms of Cartesian 
distortions is 2k x k y e xy — k y 2 e xx — k x 2 e yy = 0. In terms of 
physical strains of (2.5) and with only one shear s = 3, this 
becomes the single compatibility constraint as in (2.7), 

dei + 2 e 2 + 3 e 3 = 0. (Al) 

For the square lattice, the compatibility coefficients are 

Oi(fc) = -j^k 2 ,0 2 (k) = +^(k x 2 - k y 2 ),0 3 (k) = 

2k x k y , while for the triangular case, 0\{k) = —k 2 ,0 2 (k) = 
(k x 2 — k y 2 ),0 3 (k) = 2k x k y . In the displacement represen- 
tation, for a 2D square lattice there are two independent vari- 
ables (u x ,u y ) per lattice point (or per unit cell). In the equiv- 
alent (symmetric) strain representation, there are three phys- 
ical strains ei,e 2 ,e 3 and one constraint, so there are again 
3 — 1 = 2 independent variables. We now derive compatibil- 
ity kernels for transitions with increasing numbers of non-OP 
strains, n = 1, 2. 

1. n = 1 cases: 

Triangle/center rectangle, square/oblique, and trian- 
gle/oblique transitions: 

For these three transitions, the two-component OP strain is 
(s 2 ,e 3 ). The single non-OP strain is e\, and the harmonic 
term is f non = Ai\ei\ 2 . Substituting from compatibility 
e\ = — J2e= 2 3 Oiei jO\ , immediately yields the 2x2 matrix 
kernel of components' 3 , with prefactor v(k) inserted, namely 
Uw = vOiOtIO 2 so 

U 22 = v0 2 2 jO x 2 , U 33 = v0 3 2 IO x 2 
U 23 = v0 2 3 IO x 2 = U 32 . (A2a) 

Or explicitly, 

U 22 = v{k 2 - k y 2 ) 2 /k\ U 33 = iy{2k x k y ) 2 /k 4 

U 23 = v2k x k y (k 2 - k y 2 )/k 4 = U 32 . (A2b) 

Although these 2D transitions all have Nop = 2 and 
the same compatibility kernel, they of course differ in their 
Landau or coefficients, and in the different Ny values of 
their nonzero pseudospin vectors (pointing to corners of a 
triangle, square, and hexagon, respectively). 

2. n = 2 cases: 

(a) Square/rectangle transition: 
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The single OP is the deviatoric strain e 2 . The non- 
OP strains are the 2D compression and shear e-y and 
e 3 , and the harmonic non-OP free energy is f non = 
^ i=1 3 (Ai/2)|ei(fc)| 2 . The compatibility condition yields 
e 3 = -Ea=i2^a e <" where O a = O a /0 3 . Substitut- 
ing yields f non = (Ai/2)[|ei(fc)| 2 + J2 a ,0=i, 2 G a ,pe a ep*], 
where G aJi = {A^/A^O^p. 

Freely minimizing in the remaining non-OP variable e\ as 
df n on/del(k) = 0, yields e\{k) = B 12 (k)e 2 (k). Substi- 
tuting back, the local non-OP term f non becomes a nonlocal 
compatibility potential for the OP distortions f C ompat(^2) = 
(A 1 /2)U(k)\e 2 (k)\ 2 , where the kernel is 

U(k) = {G 22 + R 22 , n )/[l + G 11 ], (A3) 

where i?22,ii = G 22 G\\ — G\ 2 2 ■ The structure is similar to 
the 3D cases below. However, for the 2D case, the remainder 
term i?22,n = so this becomes 

A 1 U{k) = ^0 2 2 /[{Oi 2 Mi} + {0 3 2 /A 3 }]. (Ma) 
Or explicitly, 

ArUik) = lyA^kJ-ky^/i^ + iSA^A^ky) 2 }. (A4b) 

We fix 2A\jA 3 ~ 1 for simplicity, so the strength of the 
compatibility potential is determined by A\, that is essen- 
tially the elastic anisotropy parameter^ by (2.18d). Plots of 
U have been given elsewhereJiia. The same result is ob- 
tained through Lagrange multipliers: minimizing {/ non — 
A (E a =i,2,3°a e a)} in e i, e 3 yields e l = AOi/Ai for i = 
1,3. Demanding compatibility fixes the Lagrange multiplier 
A = -0 2 e 2 /[{Ol/Ai} + {Ol/A 3 }], so ei, e 3 are in terms of 
e 2 , yielding the same kernel (A4). 

As mentioned in the text, the square/rhombus is not distinct 
from the square/rectangle transition. Nonetheless, the OP is 
now e 3 and the non-OP are ei,e 2 . The kernel is just a2n3 
label interchange, 



The 3D compatibility conditions of (2.3) in Fourier space 
for Cartesian distortions Jnc(e) = k x e(fc) x k = 0, can 
be written as three equations from diagonal components of 
Incompatibility, like Inc(e) xx = 0, 

2kyk Z €y Z ky€ ZZ -\~ k z Gyy, 

'2k z k x e zx — k z c xx -t- k x c zz , 

"2k x kyG X y — k x Cyy ~\- kyC xx . (-AT) 

There are also three equations from the off-diagonal compo- 
nents, like Inc(e) yz = 0, 

kyk z c xx — k x €y Z -\- k x kyC zx -t- k z k x c X y , 

k Z k X Cyy kyC ZX -f- kyk Z C X y ~j~ k X ky&y Z , 

k x ky€ zz kzc x y -f- k z k x Cy Z -\- kyk z c zx , (A8) 

but these are not new constraints: solving (A7) for shears and 
substituting, (A8) is satisfied as an identity. In the displace- 
ment representation, for a 3D cubic lattice there are three 
independent variables (u x ,u y ,u z ) per lattice point (or per 
unit cell). In the equivalent (symmetric) strain representation, 
there are six physical strains e±, .., eg and three constraints, so 
there are again 6 — 3 = 3 independent variables. 

The 3D St Venant constraints of (A7) in terms of the Carte- 
sian distortions can be written in terms of the cubic-lattice 
physical distortions (2.6a), by inverting the coefficient matrix, 
to get 



A-JJ$) = uO z 2 /({O l 2 /A 1 } + {0 2 2 /A 2 }). (A5) 
(b) Rectangle/oblique: 

The OP is again e 3 but the non-OP energy f non = 
^ ± (A±/2)|e±| 2 is now harmonic in the combinations e + = 
+ e 2 ) and e_ = i(ei — e 2 ). The substitution for e\ 
through compatibility, and free minimization in e 2 yields the 
kernel with 0± = 0\ ± 2 as 

AtU(k) = V O s 2 /[{0 + 2 /A+} + {0- 2 /A„}]. (A6) 
B. 3D transitions 



eyy = Tf 1 " 7f 2 + 7f 35 

e„ = - -^e,. (A9) 

This yields the 3D St Venant constraints in terms of physical 
distortions, 

0[ s) ei + { 2 s) e 2 + 0i 3) e 3 + { s s) e s = 0, (A10) 
labelled by the three shears, s = 4,5, 6. Defining 



We use the 3D compatibility constraints to obtain St Venant 
kernels for n = 3, 4, 5. 



0^ = 0^/O s , (All) 
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the shears es, e 6 are related to the non-shears ei, e2, e 3 by 

e. = - ]T Oi s) e Q . (A12) 



a = l,2,3 



The cubic-lattice compatibility coefficients O a of (A10) 
are evaluated from (A7) and (A9) as 



O 



d) 



-1 



2 , i.2\ ^,(4) 



V3 



^(2fc2 - fcj), of = Oi = k y k z . (A13a) 
v6 



V2 



O 



(6) 



o, o, 



(6) 



6 = fc z V (A15c) 



These will be used for the tetragonal/orthorhombic transition, 
below. 

It is useful to define a dimensionless variable, analogous to 
the 2D version above, namely 



G aP = £ (A s /Ax)Of Oj, s) ,(A16a) 

s =4,5,6 

where G Q( g = Gp a is symmetric, and a remainder term as 
before is defined through products, 



R a p,jS = G a pG 1 s — G ai Gj3s- (Al6b) 



of = —{kl + kl), o ( 2 



2 , 1,2 N n (5) 



V3 



V2 



{ £ ] = -±=(2k 2 x -k% of = 5 = k z k x . (A13b) 
v6 



of = + fcjl), O^ = 6 = fc,fc,. (413c) 

These will be used in the cubic/tetragonal, cu- 
bic/orthorhombic, and cubic/trigonal transitions, below. 

For the tetragonal-lattice physical distortions of (2.6b), the 
Cartesian components can be written as 

e xx = -4=(ei + e 2 ); e yy = —f={e\ - e 2 ); e zz = e 3 , (414) 



,2 , 7.2N ^,(6) _ 



V2 X 



and the tetragonal-lattice compatibility coefficients from (A7) 
and (A 14) are 



1 "Tl 7 ' 2 "71' 



Then the kernels for the four 3D transitions are obtained 
in terms of the G a p, by a substitution/minimization method 
similar to 2D. 



3. n = 3 case: 

Cubic/trigonal transition: 

For this transition, the OP strains are the shears {e^} = 
{e s } = e<4, es, e§. The non-OP strains are e\, e%, e^, and the 
harmonic term is f non = J2i=i 2 3(A/2)|ei(fc)| 2 , with de- 
viatoric coefficients A2 = A3 by symmetry^. Compatibility 
here gives the shear OP e^, es, e 6 in terms of the non-OP as in 



(A12), 



E 



i=l,2,3 



Oj e, or in matrix form, 




(4) 



Of Of 



(5) q(5) ^(5) 
(6) 



0} 

V^ 6) of o 3 




(417) 



The non-OP can be written in terms of the OP shears by in- 
verting the 3x3 coefficient-matrix M above to directly yield 
^ = V. />',.' ,. Here B = M~ x = (2/DetM)N, 
where N = adj(M)/2. The determinant can be written as 

DetM/2 = [Of N 14 + 0[ 5) N 15 + Of N le \. The compo- 
nents Ni s can be evaluated, such as 



Of = -k 2 y , of = 4 = k y k z . (Alba) 



Of = of - zM 



V2 



o 



(5) 



fc2, of = 5 = fefe. (4156) 



>/2' 



01 



(G) 



(6) 



1 



-(kl-kl-kl). (418) 



2V304 

The elements A^4, iVis, -^16 j •■•A r 36 in matrix form are: 



2\/304 
(fc|+4) 
2V2O4 
{ 2kl+kl-kl 

\ 2Veo 4 



_ 2\/205 
(2fc 2_ fc 2 +fc 2 

2Veo 5 



2\^30e 2 

2V2O6 
(2kl + kl+k 2 y ) 

2%/60 6 



(^19) 
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Evaluation of DetM with (A19) yields the simple result 
2/DetM = +1, so B is = N is , and so the non-OP in terms 
of the OP strains are = 2~Z S =4 5 6 Ni S e s . 

The St Venant compatibility terms are ob- 
tained by simple substitution into the non-OP har- 
monic terms, f non = £4=1,2,3 (-"V 2 ) I e »(fc)| 2 = 
Y,e t e = 4 >5 ,6(Ai/2)Uu>(k)ee(k)e e >(k)*, where the 3 x 3 
compatibility matrix kernel for the cubic/trigonal transition is 

A l U u >(k) = v A i N ie N w, (-420) 

i=l,2,3 

and can be numerically evaluated in simulations. 
4. n = 4 case: 

Cubic/ tetragonal and cubic/ orthorhombic transitions: 
In both cases, the two-component OP are the two 3-D de- 
viatoric strains (e 3 , e 2 ). The non OP are the remaining com- 
pression and shear strains ei,e4,e5,e6, and their harmonic 
terms are f non = (A 1 /2)\e 1 \ 2 + E s =4,5, 6 (4 s / 2 )|e s (fc)| 2 , 
with shear coefficients A4 = A$ = A§ by symmetr y 14 ' 16 . 
From the compatibility equations (A12), we have e s = 
— 5Zq=i 2 3 e a and hence, using the definition (A16a), 



there are three Lagrange multipliers A' s '. Then one finds 

e s = A<»0 S /4 S ; ei = £ s AWO^/^i. From com- 
patibility, A« = -(As/O s )Y,, =2 ^[0 [ ;\l + G n ) - 

o[ 8 ^Gu]ee/(l + Gn), yielding the same kernel as before. 
This kernel was used earlier~, but is here and in^ explicitly 
stated. 

5. n = 5 case: 

Tetragonal/ 'orthorhombic transition: 

For this transition, the single-component OP is one of the 
deviatoric strains, e%. The non-OP strains are ei, e 3 , e±, es, e^, 
and the harmonic term is f non = (Ai/2)|ei| 2 + (A 3 /2)|e 3 | 2 + 
£ f =4,5,6(4 s / 2 )M 2 , with A 5 = A 6 by symmetry^ 6 .. Substi- 
tuting with (A12), but now with the compatibility coefficients 
(A15) for the tetragonal case, 



Uon= -y[ (AiMi)|ei| a + G a , e a ep% (A24) 

i=l,3 cc,/9=l,2,3 

where the G a p is defined in (A16a). 

Minimizing freely in ei,e3 and inverting a 2 x 2 matrix 
yields 



fnon = {A 1 /2)[ | ei | 2 + ]T G a> pe a e *}. {A21) 

a, 0=1,2,3 

Minimizing, we get e\ = ^£=2 3 Buei with Bu = 
—Gu/(1 + Gn) similar in structure to the square/rectangle 
case. It is easy to check using (A16a) and (A13), that Gn 
depends only on cubic-invariant combinations of k x ,k y ,k z . 

Substituting back into (A21), f non (ei) = / CO m P at(e 3 , e 2 ) 
where the St Venant term is 

fcompat= ( A i/2)U ul (k)e e (k)e e (k)*. (A22) 

i,i'=2,3 

Here, the compatibility kernel for the cubic/tetragonal and the 
cubic/orthorhombic transitions, is the 2x2 matrix 

Uw = v[G u > + Ru',u]/(1 + G n ), (423) 

with Rw,u = GwGn — Gi\Gm- This kernel can be 
numerically evaluated in simulations. 

As in all cases, the same results can be obtained by 
minimizing {/„„„ - £ s =4,5,6 A(s) E a °* )e a}} where 



ei = -[(4 3 /Ai)(?i2 + i?33,i2]e 2 /Go; 

e 3 = -[Gaa + flii,3a]e 3 /Go, (425a) 
where R aa ,p-y is defined in (A16b), and 

Go = [({As/A,} + G 33 )(l + Gn) - G 2 3 ]. (A25b) 
The kernel for the tetragonal/orthorhombic transition is 

ArUik) = iy[(A 3 /A 1 )G 22 + T 2 ]/G , (426). 

With T 2 — (4 3 /^4i)i?22,ll + -^22,33 + {G22-R33,ll — 

Gi2i?33,i2 — G 3 2i?n, 3 2}- The kernel can be evaluated nu- 
merically in simulations. 

As a check, we take uniformity in the z direction, or k z —> 
0, when 04,05 —> 0, and it is clear from compatibility that 
e 3 = e4 — e 5 = 0. Then one recovers, with 4A 6 —> A3, 
precisely the form of the square/rectangle kernel of (A4). 

Finally, the compatibility kernels in Fourier space are plot- 
ted for the four 3D transitions in Fig. 4 below, reflecting the 
high temperature unit-cell symmetries. 
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FIG. 1. Schematic figures of five distinct ferroelastic transitions 
driven by strains in 2D, with lower symmetry variants on the right: 
square to rectangle (SR); rectangle to oblique (RO); triangle to 
centered-rectangle (TR); square to oblique (SO), triangle to oblique 
(TO). For the TR case, we mark one of three equivalent lattice-point 
boxes, that becomes the new centred-rectangle unit cell of the first 
variant on the right, under the simplest distortion order-parameter. 
The other two equivalent variants have similar unit-cells, that are 
simply from ±27r/3 order-parameter rotations of that distortion. For 
the TO case, the equivalent rotations are integer multiples of 2tt/Q. 



FIG. 2. Schematic figures of four ferroelastic transitions in 3D, 
with the lower symmetry variants on the right: (a) tetragonal to or- 
thorhombic; (b) cubic to tetragonal; (c) cubic to trigonal with two 
other variants not shown; (d) cubic to orthorhombic with four other 
variants not shown. 



FIG. 3. Scaled free energy versus Nop order parameter compo- 
nents, with minima at austenite zero state and Nv martensite vari- 
ants, and with arrows to minima denoting pseudospin vectors, (a) 
Minima on a line for Nop = 1, Nv = 2, square to rectangle (and 
also tetragonal to orthorhombic, rectangle to oblique); (b ) Minima 
on a triangle for Nop — 2, Nv = 3, cubic to tetragonal (and tri- 
angle to centered rectangle); (c) Minima on a square for Nop = 2, 
Nv = 4, square to oblique polygon; (d) Minima on hexagon for 
Nop = 2, Nv = 6, for cubic to orthorhombic (and triangle to 
oblique); (e) Minima (schematic) on a tetrahedron for Nop = 3, 
Nv = 8, cubic to trigonal. 
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FIG. 4. Compatibility kernel components Ui ti (k x , k y ,k z ) in color 
plots versus (k x ,k y ,k z ) for 3D transitions, (a) Tetragonal to or- 
thorhombic case with U (k). The strength of the kernel is represented 
by a color coding in which dark brown represents the relative posi- 
tive maxima, and dark blue the minima with zero values.The z axis 
is vertical and the projections shown are at k x — 0, k y = and 
k, — 0. The maxima shown appear in the plane k y = 0, and minima 
for k x = and k z = 0. (b) Cubic to tetragonal (and also cubic to 
orthorhombic case) kernels, (i) U22(k), (ii) Us3(k), and (iii) U2s(k). 
The three components acquire positive values, and assume clover- 
leaf anisotropy in the 2D planes k x — 0, k y = and k z — 0. (c) 
Cubic to trigonal case, with a kernel component Um(k). In addition 
to the clover-leaf pattern in the plane k y = 0, there is a butterfly 
anisotropy in the plane k x — 0, similar to that of the square to rect- 
angle kernel in 2D. 



TABLE I. Generic numbers for the scaled Landau free energies, 
listed in order of increasing number of structural variants Nv- The 
columns are: (1) the type of transition; (2) spatial dimensionality; 
(3) order-parameter (OP) dimensionality Nop or number of pseu- 
dospin components ; (4) number of free energy minima at transition 
or number of pseudospin- vector states, Nv + 1; (5) maximum-order 
invariant p max retained in the free energy; (6) number of material co- 
efficients Nmat, and type of scaled quasi-universality (q-u) as in text; 

(7) condition obeyed by OP-space angular location 4>m of minima ; 

(8) 'polyhedron' from minima in Nop dimensions. 
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